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Abstract 

The pomeron model of Donnachie and Landshoff does not conserve the electromag- 
netic current when applied to diffractive reactions such as electroproduction of a quark- 
antiquark pair or of a vector meson. We propose a treatment of this problem which 



£NJ ■ ensures a physical behaviour of cross sections in the photoproduction limit and show 

that it leads to results rather similar to those obtained from two-gluon exchange. 

OO 

*q} 1 Introduction 

The description of the pomeron within the framework of QCD remains one of the great 
tasks in strong interaction physics. Important progress has been made in this field over the 
last years, especially in the study of diffractive processes with a hard scale, for instance a 
large photon virtuality Q 2 in ep collisions. On the other hand, phenomenological models are 
still of interest in this new dynamical regime: calculations in the framework of multi-gluon 
exchange can be of considerable complexity, and it is useful to have models that reproduce 
their results in a simpler way, thus allowing the investigation of more complicated reactions 
with a reasonable amount of effort. Furthermore, such models provide an opportunity to 
extrapolate from the hard to the soft diffractive regime, e.g. to take the limit of small Q 2 in 
diffractive ep scattering, and to make contact with elastic and diffractive hadron scattering, 
domains where the application of QCD is much more difficult since perturbation theory is 
of little help. It goes by itself that one cannot expect such models to reproduce all results 
and aspects of more sophisticated and more fundamental QCD approaches. 

In this paper we are concerned with the pomeron model of Donnachie and Landshoff 
(DL) m and its confrontation with the QCD motivated model of nonperturbative gluon 
exchange by Landshoff and Nachtmann (LN) @, |3[] . It was found in that the DL model 
has problems with electromagnetic gauge invariance, which appear for instance when it is 
applied to the production of a quark- ant iquark pair in diffractive ep scattering. In Sect. |2| 
we explain the origin of this problem and propose a possible solution which makes the DL 
model reproduce several, although not all results of the two-gluon calculation in the LN 
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approach. We then investigate exclusive vector meson production under the same aspects. 
In Sect. [| we discuss our way to handle the gauge invariance problem from the point of 
view of contact terms, and make some comparison with other approaches in the literature. 
We summarise our results in Sect. ||. 

2 Diffractive (/(/-production 

2.1 The problem in the DL model and a solution 

In the DL model the pomeron couples to quarks via their vector current, a choice originally 
motivated by the analysis of elastic and diffractive hadron-hadron scattering ||. In this 
sense the pomeron is said to behave like a photon, except that it has charge conjugation 
parity C = +1 unlike the photon with C — — 1. To implement this difference of quantum 
numbers the model prescribes to subtract Feynman diagrams which are related by reversing 
the charge flow of quark lines, instead of adding them as one would do if the pomeron were 
replaced by a photon. 

Calculating the forward 7*p scattering amplitude at high energy in this model DL found 
that in order to obtain Bjorken scaling of the proton structure function at small x the quark- 
pomeron vertex cannot be pointlike but must be softened at large quark virtualities 
They made the ansatz 

PoYf(kl-m 2 q )f(k 2 2 -ml) (1) 

with 

f(k 2 ) = ~ k ") (2) 

for the pomeron coupling to quarks of mass m q and momenta k\ and k 2 , with the constants 
(3o ~ 2 GeV -1 and /xq ~ 1 GeV determined from phenomenology. 
We now investigate diffractive gg-production, 

7*(?) +P(P) - q(P q )+q(P 9 ) +P(p') , (3) 

with four-momenta given in parentheses. We will further use the variables A = p — p', 
t = A 2 , W 2 = (p + q) 2 , Q 2 = -q 2 , M 2 = (P q + P q ) 2 , (3 = Q 2 /(2A-q) and £ = (A-q)/(p-q). 
In the DL model this reaction is described by the two diagrams of Fig. |T], taken with 
opposite relative sign. One finds |4j that the electromagnetic current is not conserved here: 
the 7*p cross section ajj for the unphysical photon polarisation £3 = j Q does not vanish 
but instead has a 1/Q 2 singularity in the photoproduction limit. In a technical sense 
the reason is that the two diagrams in Fig. [1| have the "wrong" relative sign which was 
assigned in order to implement the correct charge conjugation parity of the pomeron — for 
photon instead of pomeron exchange gauge invariance is of course guaranteed. If one uses 
the Feynman gauge for the photon field and defines longitudinal and transverse photon 
polarisations with respect to the 7*p axis in the cm. of the collision then the cross section 
a i, for longitudinal photons has the same unphysical 1/Q 2 behaviour at small Q 2 , only the 
transverse cross section a? is not affected by this problem. 
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Figure 1: Diffractive gg-production in the DL model. The relative sign of the diagrams is 
the opposite of what it would be if the pomeron were replaced by a photon, because the 
pomeron and photon have opposite charge conjugation parity. 

It is instructive to look at the dependence of the cross sections on the transverse mo- 
mentum p t of the produced quark in the 7*p cm. and on the invariant mass M of the 
gg-pair. Taking quark masses equal to zero the unphysical and longitudinal cross sections 
have a factor p 2 t /M 2 relative to the transverse one, and as a consequence one finds that 
at large Q 2 the p r integrated cross sections da v j \dM 2 ) and da^j \dM 2 ) are suppressed by 
f^o/Q 2 ' l°g(Q 2 1 Vo) compared with daxj [dM 2 \ In this sense one can say that the gauge 
violating terms are of higher twist and that the model can be used for leading twist quan- 
tities. In jet production, where p 2 /M 2 is not small, the unphysical cross section is however 
not negligible compared with dar/ (dp 2 dM 2 ), even at large Q 2 , and the model is of no use 
as it stands. 

Let us have a closer look at the origin of the bad Q 2 -behaviour of dai/(dp 2 dM 2 ). From 
the expression of the longitudinal polarisation vector in Feynman gauge 



it is clear that the part of e proportional to p gives a factor Q in the amplitude and therefore 
a factor Q 2 in the 7*p cross section as required by gauge invariance. The unphysical 
behaviour comes from the part proportional to g, which would give no contribution if the 
electromagnetic current were conserved in the model. 

Without current conservation the results of a calculation will clearly be gauge depen- 
dent. Therefore we need to fix a photon gauge in order to make the model well defined. We 
feel that this is legitimate in the context of a phenomenological model, a choice of gauge 
then has to be justified by its results. We choose to work in a noncovariant gauge with 
gauge fixing vector p, where the photon field satisfies A -p — 0. In this gauge a polarisation 
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vector e in Feynman gauge becomes 



— Q • (5) 
p-q 



It is instructive to compare the tensor structure of the photon propagator, 

gT = -sgn(g 2 ) efc£ - e?e£ - + sgn(g 2 ) e^ 3 (6) 
in Feynman gauge and 

ST ~ — W + flV ) + 7^ 9*V = -sgn(g 2 ) effeg - e^i - e£e£ (7) 

in our noncovariant gauge, where £i and £2 are two orthogonal vectors transverse to p and 
q. We see that in the amplitude for the electroproduction process ep — > ep + qq the term 
£363 on the r.h.s. of drops out when the index \x is contracted with the electron current 
since this current is conserved, so that we are left with the ill behaved contribution from 
the contraction of q v /Q in Eq with the current of the produced quarks. In the gauge (|) 
the transverse polarisations, which dominate in the p f -integrated cross sections at large 
Q 2 , are the same as before, S\ = E\ and £2 — £2, whereas the longitudinal polarisation 

e = [p q] (8) 

s/l+P 2 Q 2 /(p-q) 2 P-Q V P'Q J 

behaves like Q in the photoproduction limit and leads to a reasonable behaviour of the 
amplitude. From the l.h.s. of (0) we see that when changing from Feynman to our nonco- 
variant gauge we have effectively added the contraction of —p^q v j (p ■ q) with the leptonic 
and hadronic currents to the ep — > ep + qq amplitude, the terms with g M in the photon 
propagator giving zero when contracted with the electron current. 

Instead of (|5]) one may consider other noncovariant gauges A ■ n = with a suitably 
chosen vector n. Introducing longitudinal and transverse photon polarisations with respect 
to n and q and repeating the arguments that lead from (|4]) to @ one will find that the 
ep amplitude has a well behaved photoproduction limit. This amplitude will in general be 
different for different choices of n since one has effectively added —n^q u / (n-q) to the tensor 
g^ u in the photon propagator. Note that the 7*p amplitude for photons that are transverse 
with respect to p and q will in general also be modified by such a choice of gauge. 

A candidate gauge fixing vector is n = A, which for t = is equivalent to n = p since 
then A = £p, but leads to different results at finite t. Taking n = A looks more symmetric 
with respect to the photon dissociation and proton scattering parts of our reaction; to 
leading order in £ -1 the choice n = p is however equivalent to n = p + q, which is also 
symmetric. We shall come back to the question of p versus A in Sect. ^. 

2.2 Comparison of the DL and LN models 

We shall now see that in the gauge (0) the predictions of the DL model for our process 
are remarkably similar to those in the LN model, although they are not identical. We 
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restrict ourselves to t = but admit a finite quark mass m q in the produced gg-pair. With 
Hand's convention for the flux factor the cross section of the process (§) for a given photon 
polarisation reads 



da 



dtdpt dM 2 



e 



t=0 



1287T 3 (Q2 + M 2)2 M 2 1 1 _ 4 (p2 + m 2) / M 2 



e)\ 2 



(9) 



spins 



where J2's P ins stands for spin summation for the final state particles and spin average for 
the initial proton. In the DL model one has 

16iree q j-q 



M(e) 



r Qjp(0) ■Ku{P q ) 1 -ev{P q ) 



with 



K 



3 2p-q 

9(3 2 f{i-m 2 g ) + f(u-m 2 q ) 



(10) 



'11^ 



4vr 2 

Here t = (q — P q ) 2 and u = (q — P q ) 2 are the Mandelstam variables of the pomeron-photon 
subreaction, e q is the charge of the produced quark in units of the positron charge e and 
ajp(t) is the soft pomeron trajectory. 
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u{p') 



F 1 {t) 1 »--—a^KF 2 {t) 
2m n 



u(p) 



(12) 



is the isoscalar nucleon vector current with F\ and F<i being the isoscalar Dirac and Pauli 
nucleon form factors, i.e. the sum of the respective form factors of proton and neutron. 

In the LN model the reaction (^) is described by the diagrams of Fig. |2| plus those 
obtained by reversing the charge flow of the upper quark line; for details of the calculation 
cf. PJ. Key quantities of the model are the moments 



dl 2 \af D(-l 2 )} 2 
dl 2 [a^D(-l 2 )] 2 -l 2 



9/3, 



47T 

8tt 



(13) 



of the nonperturbative gluon propagator D(l 2 ), where the strong coupling in the nonper- 
turbative region is taken as ~ 1 following ||. (3$ and /io have been identified in || 
with the corresponding parameters in the DL model, cf. (|l|) and (0). Note that /Zq has here 
the significance of a typical gluon virtuality I 2 dominating in the integrals. For diffractive 
gg-production we have in our noncovariant gauge 



M(e) 



lQ7iee q j.q 



2p-q 



-ajp(0) 



\ 



a s (X 2 

(0) 



L 1 u{P q ) 1 -ev{P q ) 



- (Li - L 2 



A-q 
1 



u{P q ) 1 -Av{P q 



+ (Li - L 2 ) — m q u(P q ) 1 -e 1 -Av(P q ) 



(14) 
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where _ 

poo 

Li= dfia^Di-l^Ml^plX 2 ) , * = 1,2 (15) 

JO 

are loop integrals with weighting functions fi, f 2 whose expressions can be found in [|J. 
The variable 

A-4±^ (16) 

can be seen as the relevant hard scale of the process (7|, [|. The square root of a s (A 2 )/a^ 
in ( |14|) corresponds to taking the quark-gluon coupling at a perturbative scale only for the 
upper left vertex in the diagrams of Fig. where one quark leg has a virtuality of order 
A 2 . This choice was made in |1 and [||, JTUL p~T|. One may argue that the coupling should 



be taken as perturbative at both upper vertices, then there is no square root in (0), (|I7|). 
A similar comment holds for vector meson production (PU|). 

At this point we remark that many features of the LN model are due to the two-gluon 
exchange mechanism and in common with other two-gluon approaches. The authors of H 



and |T2| use the gluon density in the proton to describe its coupling to the two exchanged 



gluons; at t — their expressions are related to those in the LN model by the substitution?] 

71 d 



A Qs ' dP 



\ 



a,(A 2 ) . £1-^(0) [a (o) D( _ Z 2 )]2 . Z 2 f (17) 



(0) 

as 



where <?(£, I 2 ) is the gluon density at a factorisation scale I 2 and a s on the l.h.s. taken 
at scale A 2 in || and I 2 in [12|]. The principal difference in the predictions of the two 
approaches is thus the ^-dependence and overall normalisation, whereas the dependence 
on Q 2 , M 2 , p\ and the quark mass comes out very similar. 

Coming back to the calculation in the LN model one finds that the contribution of 
diagram (a) in Fig. |2] plus the one with reversed quark charge flow is obtained by replacing 
both Li and L 2 with 9/3q/(47t), so that only the first term in the brackets of (|14"|) survives. 
Except for the root of a s (\ 2 )/a^ this is just the DL model expression (PUf), (|TT]) with 
the form factor (fj) set equal to 1. In the DL model this form factor ensures the decrease 
of the amplitude at large p 2 , while in the LN model the same effect is achieved by the 
contribution of diagram (6) and its analogue with opposite quark charge flow. In the LN 
model one can of course also calculate with the polarisations e in Feynman gauge instead 
of e since the sum of the four diagrams is gauge invariant. If one does this M then the 
diagrams of type (a) give contributions the amplitude for longitudinal photons that do not 
vanish at Q 2 — > 0, which are exactly cancelled by the diagrams of type (b). In the DL 
model the global factor [f(t — m 2 ) + f(u — m 2 )]/2 cannot achieve this and one is left with 
an unphysical behaviour at Q 2 — > when working in Feynman gauge. 

Comparing (|T0|) with (|14]) we see that the expression in the DL model can be obtained 
from the LN results by replacing both loop integrals L\ and L 2 by K and by dropping 
the ratio a s {\ 2 ) / af> of strong couplings at different scales in the LN expression. This 

1 This substitution is consistent with M, where the gluon density was estimated within the LN model. 
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Figure 2: Two of the four diagrams for gg-production in the LN model, the other two are 
obtained by reversing the charge flow of the upper quark line. The lower line stands for a 
constituent quark in the proton and the pomeron is approximated by two nonperturbative 
gluons, represented by dashed lines. 



ratio gives a different normalisation as a function of the scale A 2 , whereas the differences 
between the loop integrals ( f[5|) and the combination of form factors ( |TTD have more subtle 
effects as we shall now see. 

At t = the longitudinal polarisation vector reads eft = Q/(A- q) ■ A M to leading order 
in and it is easy to see that the sum of the first and the second term in the brackets 
of ( |I3D gives an amplitude proportional to L 2 for longitudinal photons. For transverse 
photons the amplitude is proportional to L\ in the case m q = when only the first term 
contributes, and it involves both loop integrals for finite quark mass due to the third term. 
These differences are all absent in the DL result, where the effect of the form factor is the 
same for all photon polarisations. 

To illustrate the differences let us compare the values of L 1; L 2 and K in some kine- 
matical limits, keeping in mind that there is also an additional root of a s (\ 2 )/a^ in (|T4]). 
We always assume that /j% <C Q 2 + M 2 so that from (fj) , ( |TT| ) °ne has 



K _ 9$ t4 



8tt /i 2 + A 2 

up to terms suppressed by Hq/{Q 2 + M 2 ). 

For A 2 ^> /i 2 , a reasonable approximation of the loop integrals L\, L 2 is obtained by 
expanding the weighting functions f±, f 2 around I 2 = 0, and with (|T3| ) one finds 

^^A 2 " V AV ' 87TA 2 " \ ~ ~~>? J ' 87TA 2 " (19) 

In the case p 2 ^> m 2 , which is relevant for jet production, one has A 2 ~ p 2 /{l — 0) so 
that due to the particular choice of form factor in (EI) the p 2 -dependence of the amplitude 
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is the same in the two models for all photon polarisations. The /^-dependence is however 
different: in the LN model one finds a zero of the longitudinal amplitude around (3 = 1/2 
and a suppression of the transverse one at small /3, both phenomenologically important 



effects [pT0|, §[] which are not found in the DL model. 

In the opposite limit, p 2 <C m 2 at A 2 fi^, which is relevant for heavy flavour produc- 
tion, one has p 2 /X 2 <C 1 and thus L 2 ~ K. In the various 7*p cross sections (eq. (27) of 
[fT~0|1 ) one finds that terms going with Li are suppressed by powers of Pt/rn q compared with 
those only depending on L 2 , so that in this limit the two models lead to the same results, 
apart from the ratio of strong coupling constants. 

When both p t and m q are small so that A 2 <C /i 2 , the weighting functions f\ and f 2 tend 
to 1 and we have 

LxM^-f , ^ir> (2°) 

47T 87T 

so that in this limit the DL amplitude is just half of that in the LN model, given that 
a s (\ 2 )/a^ in the LN expression should be replaced with 1 at small A 2 . 



2.3 The longitudinal diffractive structure function in the DL 
model 

As an application of our modification of the DL model we calculate the integrated cross 
section for diffractive gg-production, which can be taken as an approximation for the 
inclusive diffractive cross section — except in the region of small j3 where diffractive states 
with additional gluons are known to be important. Expressing our result in terms of the 
diffractive structure functions and F^ 4 \ conventionally defined by 



da, 



ep 



dx dQ 2 dt; dt 



xQ 4 



l-y + y 2 /2) F T D(4) + (1 



F 2 D{4) = F^ (4) + Ff (4) 



(21) 



we obtain for Q 2 + M 2 > $ 



f F (t,t) ■ F£ L (P,Q>,t) 



(22) 



with 



and 



9ft 
4vr 2 



l-2a F (t) 



F^t) 2 - —F 2 {t) 



Am 2 



(23) 



Ff(J3,Q 2 




+ 




(24) 
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where we have summed over massless quarks u, u, d, d, s, s. We have made the approxima- 
tion t = in and but kept the strong i-dependence from the pomeron trajectory 
and the elastic form factors in fjp. The terms of higher order in ( p4"|) give only small 
corrections: we find that the effect of the terms of order P^q/Q 2 in F^f is below 3% for 
Q 2 /(3 > 5 GeV 2 and that the approximation for Ff is better than 4% for Q 2 / (3 > 10 GeV 2 . 

The result for F^ is what has long been obtained by DL |1] . It shows a scaling behaviour 
while the longitudinal structure function is power suppressed in 1/Q 2 , with Q 2 ■ Fj^ having 
only a weak logarithmic dependence on Q 2 . As (3 — > 1 the transverse structure function 
vanishes like 1 — (3 whereas the longitudinal one goes to a constant. All these features 
are also found in the LN model ||11|| . In the DL model the ratio F^ /F^ turns out to be 



numerically rather big even at moderately large (3: for Q 2 between 10 GeV 2 and 30 GeV 2 
it becomes equal to 1 at (3 between 0.65 and 0.8. Taking Q 2 = 25 GeV 2 one finds F^ 
rising all the way up to (3 = 1, with the increase of F^ overcompensating the decrease of 
Frf for (3 > 1/2. This looks difficult to reconcile with what is seen in the HERA data |T3] . 

In the LN model one also finds that F^ becomes numerically important at large (3 
where FyT goes to zero [II], but there the effect is much less dramatic. This is because 
of the extra factor of the running strong coupling at the scale X 2 = p 2 /(l — (3). It leads 
to a suppression at large (3 which is stronger for Ff than for F^ since the former is 
less dominated by small p 2 than the latter. It thus seems that the scale dependence of 
the strong coupling which appears in the two-gluon calculation leads to a more realistic 
prediction for the longitudinal structure function. 



3 Exclusive production of vector mesons 

The first process for which the DL and LN models have been compared and where their close 
relation was observed [§] was elastic electroproduction of a vector meson V = p, <p, J/ip, . . ., 

l*(q)+p(p)^V(q r )+p(p') . (25) 

It turns out that the DL model in its original form also violates electromagnetic current 
conservation in this process, as was already remarked in [13]. For the 7*p cross section one 
obtains 

^= 1 V'lAif (26) 

with 

M = «^ P. V DL (s, O • ^ , M . • (27) 



M 2p-q^ UhK ' ' Q 2 + M 2 -t + 2/j 



where M denotes the meson mass, T e + e - its decay width into e + e , and e and e' the 
respective polarisation vectors of photon and meson. The polarisation dependence is given 
by 

« UlA (P'g) (A-eQ-(A- e ) (p-e>) + (e-e') jp-q) 

F D L\ e ^ e ) = > l 2 8) 

p-q 
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where we have used e' ■ q' = but not made any assumption about e so that this expression 
holds in any photon gauge. As announced V DL does not vanish for the unphysical polari- 
sation e = £3. Working in the collision cm. and defining polarisation vectors with respect 
to the p-q axis for the photon and with respect to the p'-q' axis for the meson we find in 
particular that for the transition of a longitudinal photon to a longitudinal vector meson 
one has 

n 2 - M 2 + t 

V nT (L,L) = ^ ^—^ (29) 

DLK ' ; 2MQ K ' 

in Feynman gauge, with an unphysical behaviour at small Q. In the LN model one finds^j 



Q ' (A2) V LN <s,s> 



(0) 



48M 2 37rr e+e - j.q ffl 
e V M 2p.q^ \ 

? / A [off)]'* ((1, - A t / 2) 2 ) a ((1, + A./2)") Q2 + J[~ t ^_ Am , (30) 

where At is the transverse part of A with respect to p and q, and 

p , c n _ ip-e) (A-Q - (A-e) (p-eQ + (e-eQ (p-g) + g (p-e) (p-Q (31) 

Again this expression is valid in any photon gauge. If we require Q 2 + M 2 ^> //q and 
remember that the typical gluon virtualities in the integrals flT3| ) are of order /j,q then the 
terms after V DL and V LN in (p7|) and (p0|) are equal at t = 0; their difference will only 
become important when —t is of order /Iq. Apart from this and from the root of a s (X 2 )/af\ 
where in analogy to ( |16D one would now choose A 2 = (Q 2 + M 2 )/4, the amplitudes in the 
two models then differ only by their polarisation factors V DL and V LN . The latter has an 
extra term £ (p • e) (p ■ e') in the numerator, which was overlooked in || and recently 
reported in |TJ]. It precisely restores gauge invariance and guarantees a reasonable small-Q 2 
behaviour 

V LN (L,L) = ^ (32) 



for the transition from a longitudinal photon to a longitudinal meson. ( PU[ ) and ( |3"2"D do 
not even agree in the limit Q 2 3> M 2 , where V LN (L,L) is twice as large as V DL (L, L). 
For transverse photons on the other hand V DL gives the same as V LN . Note however that, 
unlike in the case of the diffractive structure function discussed in Sect. |2.3| , transverse 
photon polarisation in vector meson production is suppressed at large Q 2 ; with one 



finds o"l/ctt = Q 2 /M 2 for the ratio of longitudinal and transverse cross sections at all Q 2 
and t. 

If we work in the gauge A-p = then the extra term in V LN vanishes identically so that 
V DL agrees with V LN for all physical photon polarisations. Defining the DL model in this 
gauge we thus find that for Q 2 + M 2 3> /i 2 , and —t /Xq it gives the same result as the LN 

2 Our expression for M differs from those in || and || [l5| by numerical factors. Making the replacement 
( |l7| ) we agree with the result of Q . 
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approach up to the scale dependent ratio a s {\ 2 ) / af^ . This corresponds to what we found 
for gg-production in the limit m 2 q 3> p\ and A 2 ^> /i 2 , in Sect. |2]^; note that following Q 
in the calculation of meson production we have used a constituent mass of M/2 for the 
quarks and neglected their transverse momentum in the meson so that kinematically the 
two processes are equivalent. 

Let us finally compare the two gauges A ■ n = with n = p and n = A at finite t. With 
A as gauge fixing vector one does not reproduce the results of the LN calculation since the 
extra term in V LN does not always vanish. This happens when the photon is transversely 
polarised in the scattering plane while the meson is longitudinal. V LN is then negligible 

in the small-£ limit whereas V DL = —\f\~t~\ /M is not: the DL model defined in this gauge 
thus violates s-channel helicity conservation. Both from this point of view and in order 
to reproduce as closely as possible the two-gluon result the choice n = p therefore seems 
preferable to us. 



4 Discussion 

4.1 Contact terms 

The point of view underlying our detailed comparison of the DL and the LN models is that 
diffraction can be described in QCD by multi-gluon exchange in the colour singlet channel, 
and that the LN model is a simple implementation of this idea. The DL model, on the 
other hand, looks more like an effective theory in which gluon degrees of freedom are no 
longer explicitly present. 

It can be shown that in the high-energy limit the coupling of two t-channel gluons to a 
single quark line can be written in terms of the quark vector current 0. Thus the diagram 
in Fig. (a) is equivalent to the one in Fig. p] (a) calculated without a form factor for the 
quark-pomeron coupling, as we reported in the sequel of flT7|). Technically this is seen as 
follows: for the upper quark line in Fig. ||] (a) one has + m q )^ ~ 2(p ■ k)p" to leading 
order in energy, where k denotes the momentum of the quark line between the two quark- 
gluon vertices, while p is the dominant part of the gluon polarisation at each vertex and 
becomes the dominant "pomeron polarisation" in the DL model. One can also explicitly 
see how this gives the correct C — +1 quantum number of the exchange: reversing the 
charge flow of the quark line but keeping the flow of its momentum unchanged one obtains 
rriq)^ ~ — 2(p ■ k)^; the relative minus sign is precisely what is put in by hand 
between the two diagrams in the DL model. 

The diagram in Fig. (b) has a different topology than the diagrams of Fig. [I], and 
its contribution to the amplitude cannot be entirely rewritten in terms of an effective 
quark-pomeron coupling with a form factor. It is plausible to assume that a part of its 
contribution will have the structure of a contact interaction between the quark line, the 
pomeron and the photon; note that the generation of contact interactions is well known in 
effective field theories. Such a contact interaction can have a rather rich structure and we 
are not attempting here to construct a corresponding extension of the DL model, in which 
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gauge invariance is restored. Let us however make the hypothesis that such an extension 
can be found, where the sum of all diagrams for a process conserves the electromagnetic 
current. Terms with an unphysical l/Q-behaviour from the diagrams of Fig. [1] will then be 
cancelled by terms from the diagram with a quark-pomeron-photon contact interaction. 

Now the procedure of gauge fixing we proposed in Sect. [O] can be seen in a different 
light. In a gauge invariant theory one can of course work in any photon gauge; some 
terms in the amplitude will be "shifted" between different diagrams when the gauge is 
changed. There will be gauges for which the diagrams with and without contact term are 
separately well behaved in Q; an example is just the gauge A-n = with a suitably chosen 
vector n. Leaving out the contribution from the contact term (because we do not know 
its form) we then make an error, but this error is finite in the photoproduction limit. The 
error will be different for different choices of n; and whether there is a particularly good 
choice or even one where the contact term is completely eliminated we can of course not 
determine without knowing the specific form of the contact interaction. Nevertheless there 
are certain criteria: our discussion at the end of Sect. ^| shows for instance that there are 
"unfortunate" choices for which s-channel helicity in meson production is not conserved 
by the diagrams without contact terms; in a model reproducing the two-gluon exchange 
results such helicity violating terms will then only be cancelled when contact interaction 
diagrams are taken into account. 



4.2 Comparison with other approaches 

Several approaches to the problem of coupling the pomeron to quarks have been made in 
the literature. In [|17]] and later in [18| the pomeron was treated exactly like a photon, with 
a 7 At -coupling but without a form factor and without any relative minus signs between 
diagrams. By construction this model does not have any trouble with gauge invariance, 
but at the price of actually describing C — — 1 instead of C = +1 exchange. It can 
for instance not describe the forward j*p scattering amplitude, i.e. the inclusive proton 
structure function at small x, since the coupling of three vector particles (two 7* and the 
pomeron) to a quark loop vanishes due to Furry's theorem. For diffractive gg-production it 
gives a p r behaviour that is remarkably different from the one in the DL calculation. The 
DL model with the form factor (0), as well as the LN model, gives a transverse 7*p cross 
section dux / {dp 2 dM 2 ) that falls off approximately like 1/pf at large pt and is finite as pt 
goes to zero. In the approach of |T7|, |Ti| the decrease at large p t is much slower, and for 
massless quarks one has a l/p^-singularity at zero p t . This has important consequences 
for the description of the leading twist part of F 2 D in terms of diffractive parton densities, 
whose evolution equation then has an inhomogeneous term — just as in the case of parton 
densities in the photon — which is intimately related with this collinear singularity in the 
transverse qq cross section. Such an inhomogeneous term will not be generated in the 
DL model. We also remark that unlike the DL and LN approaches the pomeron model 
of [O, ITS] gives a leading twist contribution to the longitudinal structure function F®, 



going like /3 (1 — /3), in strict analogy with the photon structure function |L9 



The model of PD[ assumes a pointlike quark-pomeron coupling through the scalar in- 
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stead of the vector current of the quarks, i.e. it treats the pomeron as a scalar exchange. 
This is again gauge invariant, and the scalar current has the correct quantum numbers to 
model the pomeron, in particular it is C — +1. Taking the produced quarks as massless 
one obtains a transverse qq cross section going like 1/Pt at large and at small pt — the 
divergence will again lead to an inhomogeneous evolution equation for diffractive quark 
densities — and a longitudinal qq cross section that vanishes at t — 0. We note that for 
massless quarks a scalar quark-pomeron coupling flips the quark helicity, in contrast to 
models with a vector coupling which conserves it. The helicity of massless quarks is of 
course also conserved by multi-gluon exchange, at least if one considers the 7^-structure 
of the perturbative quark-gluon vertex to be relevant in this context. 

The starting point of |14| is a pomeron with a vector coupling and a sign factor to ensure 
C = +1 exchange; a scalar term accompanied with specific sign instructions is then added 
to the 7^-part of the coupling in order to restore current conservation in vector meson 
production. For a programme of investigating the general structure of the pomeron-quark 
vertex we refer to |2l] . 



An effective theory for high energy quark-quark scattering, starting from multi-gluon 



exchange, has been formulated in ||22|| . The current for C = +1 exchange has the structure 

ip^ D Ul D U2 . . .D U2n+l ip, where the odd number of covariant derivatives gives the correct 
quantum numbers of the pomeron. The Dirac matrix 7^ leads to the same quark spin 
structure of the coupling as in the DL model, while the derivatives guarantee the correct 
signs between diagrams when the quark charge flow is reversed, in a manner similar to the 
explicit example of two-gluon exchange discussed in the previous subsection. 



5 Summary 

There are several possibilities to describe the pomeron coupling to quarks in a phenomeno- 
logical model. The Donnachie-Landshoff model chooses a 7^-coupling and implements the 
C = +1 parity of the pomeron by introducing appropriate relative minus signs between 
diagrams. In order to obtain phenomenologically reasonable results the pomeron-quark 
coupling cannot be pointlike in this model. We point out that multi-gluon exchange in 
QCD motivates a form factor behaviour of the quark-pomeron coupling^ as well as helicity 
conservation for massless quarks, which is ensured by the 7 At -vertex in the DL model. 

The price to be paid for introducing minus signs "by hand" in this model is that the 
electromagnetic current is not conserved. We propose to define the model by specifying a 
noncovariant gauge A ■ n = in which a good Q 2 -behaviour is guaranteed for all physical 
photon polarisations in the limit where the photon becomes real. The choice of a gauge 
fixing vector n is not unique and has to be motivated, for instance by comparing predictions 
with more elaborate models such as two-gluon exchange. Our preferred choice here is n = p. 

3 This should at least hold for the part of the coupling that gives the leading twist contribution to the 
diffractive structure function F 2 D . In Q a possible pointlike structure of the pomeron in perturbative 
QCD was discussed at the level of nonleading twist. 
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A way to restore current conservation is the introduction of contact terms between 
quark, photon and pomeron lines, which seems natural if one thinks of the model as an 
"effective theory" of multi-gluon exchange. To construct such terms is beyond the scope 
of this work, but we note that in such a model one can make a suitable choice of photon 
gauge to minimise the contribution of contact term graphs, so that reasonable results may 
be obtained even when they are left out. Our choice of gauge then assures in particular a 
good Q 2 -behaviour of the diagrams with and without contact terms separately. 

We have applied the DL model in the gauge A ■ p = to diffractive gg-production at 
t — and compared in some detail its results with those of the LN model of two-gluon 
exchange. While their predictions differ in detail there is a strong similarity between the 
two models, in particular for the ^-dependence, which leads to the prediction that is 
of leading twist whereas is not. For vector meson production we find almost identical 
results in the two models when Q 2 or the meson mass is large and t is small, provided one 
takes p as gauge fixing vector. Taking A instead would lead to a violation of s-channel 
helicity conservation, at variance with two-gluon exchange. 
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